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ABSTRACT 


I  propose  a  novel  method  for  sensing  wavefront  errors  in  an  aberrated  solar 
telescope.  The  wavefront  sensor  consists  of  a  transparent  mask  located  at 
an  image  plane.  A  difference  image  of  the  scene  under  observation  is 
encoded  in  density  on  the  mask.  Wavefront  error  slopes  can  be  detected  in 
a  pupil  image  following  the  mask  in  the  form  of  intensity  variations. 

The  principle  of  the  method  is  described  and  results  of  one-dimensional 
simulation  calculations  are  presented.  The  technical  implementation  of 
the  method  is  discussed  as  well  as  the  effect  of  photon  noise  on  the  slope 
measurement.  It  is  demonstrated  that  the  proposed  method  is  very  sensi¬ 
tive  and  requires  only  a  few  percent  of  the  incident  light  for  a  photon-noise 
limited  null  state  measurement. 
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1.  Introduction 

Wavefront  sensors  generally  operate  on  confined  targets,  sometimes  assumed  to  be 
point  sources,  such  as  bright  stars  or  target  glints.  Consequently,  a  solar  adaptive  optic 
often  works  with  small,  high  contrast  structures  such  as  small  sunspots  or  pores.  The 
construction  of  a  wavefront  sensor  for  more  general  targets,  like,  e.  g.,  solar  granulation  is 
difficult  and  expensive  if  conventional  concepts  are  used.  For  example,  a  Hartmann-type 
wavefront  sensor  that  operates  on  granulation  presumably  requires  one  correlation  tracker 
for  each  subaperture  channel  and  might  be  costly. 

This  paper  describes  a  wavefront  sensor  that  is  based  on  a  novel  principle  of  wave- 
front  detection  which  is  particularly  well  suited  for  extended,  incoherent  sources  (von  der 
Liihe,  1987).  The  wavefront  sensor  consists  of  a  transparent  "mask”  in  an  image  plane  of 
the  telescope,  followed  by  optics  that  reimage  the  instrument’s  entrance  pupil  (fig.  1). 
Given  a  suitable  transmittance  distribution  of  the  mask  that  is  based  on  partial 
knowledge  of  the  fine  structure  present  in  the  target  region,  wavefront  deformations  can 
be  measured  in  the  pupil  image  in  form  of  fluctuations  in  intensity.  The  proposed  device 
results  in  a  sensitive  technique  to  measure  the  atmospheric  and  instrumental  aberrations 
of  a  solar  telescope. 

The  principle  of  the  method  is  described  in  the  next  section.  Numerical  simulations 
are  presented  in  section  3.  A  possible  way  to  realize  a  wavefront  sensor  for  a  solar  adap¬ 
tive  optic  based  on  the  technique  is  discussed  in  section  4  and  a  signal-to-noise  ratio 
analysis  for  photon  noise  is  done  in  section  5. 

2.  Principle  of  the  technique 

The  basic  optical  layout  of  the  method  is  shown  in  fig.  1,  The  target  structure  (sur¬ 
face  of  the  sun)  is  located  in  the  plane  labeled  O.  An  optical  system  that  has  an  entrance 
pupil  in  P  forms  an  image  of  O  in  the  focal  plane  O’,  Additional  optics  reimage  the 
entrance  pupil  plane  P  io  P',  Suppose  that  the  intensity  distribution  in  O  is  given  by 
/(x),  where  x  is  a  two-dimensional  coordinate.  For  the  sake  of  simplicity,  coordinates  x 
will  be  used  to  designate  object  points  as  well  as  image  points.  We  insert  a  transparent 
screen  (the  "mask")  into  the  focal  plane  O'.  The  mask’s  transmittance  m(x)  shall  be 
given  by: 

m(x)  =  B  +  C  I  /(x  -  A)  -  /{x  +  A)  I  (1) 

The  constants  B  and  C  are  chosen  such  as  to  limit  the  mask’s  transparency  between 
0  and  1.  The  fixed  shift  A  is  chosen  not  larger  than  the  correlation  scale  of  the  target 


P' 


Fig.i 

Optical  layout  of  the  proposed  method. 

structure,  i.  e.,  the  correlation  scale  of  granulation  (some  1.5  arc  seconds).  Fig.  2  presents 
a  10x10  arc  second  section  of  the  photosphere  (upper  section)  and  a  difference  picture 
according  to  eqn.  (l)  (lower  left  section)  where  the  shift  A  corresponds  to  0.5  arcseconds 
in  the  horizontal  direction. 

In  order  to  understand  how  the  technique  works,  let  us  assume  that  the  wavefront 
error  consists  of  a  purely  linear  term,  i.  e.,  there  is  only  a  displacement  of  the  focal  plane 
image.  We  place  a  detector  into  P'  that  measures  the  total  intensity  in  the  pupil  image 
following  the  mask.  If  there  is  no  displacement,  we  will  measure  a  certain  integral  inten¬ 
sity  \n  P'.  If  the  image  moves  in  the  A  -  direction,  bright  portions  in  the  held  will  coin¬ 
cide  with  sections  in  the  mask  that  have  higher  transmission,  so  the  overall  r'lcasured 
intensity  will  increase.  The  opposite  holds  if  the  focal  plane  image  moves  in  the  -A  - 
direction.  We  can  therefore  measure  overall  image  shift  parallel  to  A  by  measuring  the 
integral  intensity  in  P'.  A  second  device  having  a  mask  shift  A'  perpendicular  to  A 
measures  image  displacements  in  the  orthogonal  direction. 

This  principle  is  easily  extended  to  the  measurement  of  local  wavefront  tilt.  Con¬ 
sider  a  section  A  in  the  entrance  pupil  of  the  telescope  and  its  image  A'  in  P'  (fig.  1). 
The  section  A  will  form  an  image  of  the  target  structure  on  the  mask,  displaced  according 
to  the  local  wave  front  tilt  averaged  over  the  area  of  A.  The  displacement  can  be  meas¬ 
ured  as  a  change  in  intensity  averaged  over  the  area  A The  distribution  of  the 
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Fig.  2 


A  10x10  are  second  section  of  the  solar  surface  (top).  A  difference  mask  generated  by 
shifting  the  picture  0.7  are  seconds  to  the  left  and  right  and  subtracting  (lower  left). 
Same  as  before,  but  with  only  three  transmittance  levels  0,  0.5  and  1  (lower  right). 


wavefront  gradient  can  be  measured  by  using  a  detector  array  in  P',  and,  if  the  gradient 
in  the  orthogonal  direction  is  also  known,  the  wavefront  can  be  reconstructed  by  using 
standard  integration  algorithms. 

It  is  shown  in  von  der  Liihe  (1987)  that  the  intensity  distribution  P(x)  in  the  pupil 
image  P '  is  approximately  given  with: 

P(s)  «  P(0)  +  ABC  I  P(a/2)r*  8in(  ^(x)  -  -  «r/2)  )  (2) 

a 

X  is  a  two-dimensional  pupil  plane  coordinate.  F[t)  is  the  Fourier  transform  of  /(x). 
^(x)  denotes  the  wavefront  deformation  at  pupil  coordinate  x.  The  term  |  F(a/2)\  ’ 
represents  an  integral: 

IW2)|»  -  I  /’(x)l  >  xinsAx  d*  (3) 


HP  is  the  half  plane  in  Fourier  space  limited  by  a  straight  line  perpendicular  to  the 
A-direction  and  through  the  frequency  origin.  The  "effective  shear"  a  results  from  the 
product  of  the  power  spectrum  and  the  sine  term  in  eqn.  (3).  Fig.  3  represents  the  power 


spectrum,  the  sine  term  and  the  product. 


Fig.  3 

The  integrand  of  eqn.  (S).  (a)  represents  the  power  spectrum  of  a  bandlimited  target 
structure,  (b)  shows  the  sine  term  when  A  is  the  correlation  length,  (d)  if  ^  is 
appreciably  larger,  (c)  and  (e)  show  the  corresponding  products,  (c)  is  the  effective 
operator  used  to  generate  wavefront  differences.  With  (d),  no  local  wavefront 
gradients  can  be  obtained. 


The  shear  o  is  given  by  the  separation  of  the  two  extrema  in  the  product  shown  in  fig.  3c. 
Hence,  o  is  determined  by  both  the  equivalent  width  of  the  power  spectrum  which,  in 
turn,  is  related  to  the  correlation  scale  of  the  target  structure,  and  the  mask  shift  A.  As  a 
rule  of  thumb,  A  should  match  the  correlation  scale  of  the  target  structure  in  order  to 
make  the  device  most  efficient. 

■ 

If  wavefront  errors  are  small  compared  to  ir/3,  eqn.  (2)  can  be  linearized: 

P(.)  «  B*  F{0)  +  4BC  \F{a/2)\*  (  ^(.)  -  ^(s  -  a/2)  )  (4) 

This  equation  holds  when  an  adaptive  system  operates  close  to  zero  wavefront  error. 

In  the  real  world,  one  will  have  to  take  precautions  against  intensity  fluctuatuins 
caused  by,  e.  g.,  scintillation.  Those  would  be  interpreted  as  global  or  local  wavefront 
tilts,  produce  wavefront  fitting  errors  and  may  drive  a  corrective  system  out  of  its  range. 
One  possibility  to  eliminate  intensity  fluctuations  is  to  make  two  measurements  of  the 
pupil  intensity  with  two  masks  that  have  opposite  contrast.  The  transmitance  m '  of  the 
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inverted  contrast  mask  is  given  with: 

m'(x)  =  S  +  cj/C^  +  A)  -  /(x-A))  .  (5) 

and  the  pupil  intensity  distribution  caused  by  m'  is  for  small  wavefront  errors: 


P'(s)  «  /’(O)  -  ABC  I  F{al2)\ *  ( <6(i)  -  ^(s  -  a/2)  )  (6) 

The  difference  P  —  P '  is  independent  of  overall  fluctuations  in  intensity: 

P(i)  -  P'(*)  «  8PC  |P(a/2]r*  (  m  -  H*  -  )  (7) 


3.  Numerical  simulations 

The  performance  of  the  method  is  demonstrated  in  one-dimansional  simulation  cal¬ 
culations.  Fig.  4a  presents  a  one-dimensional  object  distribution,  the  pattern  was  chosen 
to  model  solar  granulation  with  8%  contrast. 

TEST  OBJECT 


Image  coordinate 

Fig.  4a 


The  sampling  matches  the  resolution  of  a  70cm  telescope  at  a  light  wavelength  of  SOOnm, 
the  field  corresponds  to  9.6  seconds  of  arc.  Fig.  4b  shows  three  sample  waveforms  used  to 
aberrate  the  focal  plane  image.  The  dashed  line  represents  a  quadratic  wavefront  error 
term  (defocus),  the  dot-dashed  line  a  third-order  coma  term  and  the  solid  line  a  random 


-0.500  -0.300  -0.100  0.100  0.300  0.500 

Pupil  coordinate 

Fig.  4b 

wavefront  aberration  with  second-order  statistics  corresponding  to  those  caused  by  atmos¬ 
pheric  turbulence.  Fig.  5  shows  point  spread  functions  for  all  three  wave  forms  and  for 
values  of  rms  wavefront  error  ranging  from  minus  to  plus  two  waves.  Fig.  6  presents  the 
results.  The  first  column  of  surface  plots  represents  the  one-dimensional  pupil  plane 
intensity  as  a  function  of  rms  wavefront  error,  for  all  three  waveforms.  The  second 
column  presents  the  intensity  difference  (eqn.  (7))  for  mask  measurements  with  inverted 
contrast.  The  third  column  presents  the  "reconstructed"  waveform,  calculated  by 

integrating  the  differences  from  the  second  column. 

It  is  apparent  that  the  signal  in  the  pupil  image  is  relatively  low,  it  is  of  the  order  of 
the  contrast  of  the  target  structure  (first  column).  The  pupil  intensity  difference  (second 
column)  reveals  the  non-linearity  of  the  response,  as  expressed  in  eqns.  (2)  and  (6).  The 
reconstructed  waveforms,  however,  still  have  much  resemblance  with  the  input 
waveforms,  even  for  large  rms  wavefront  errors,  except  for  sections  where  the  wavefront 
gradient  is  very  large.  One  should  keep  in  mind,  however,  that  the  "reconstructed"  wave- 
front  would  be  the  initial  response  of  an  adaptive  system  to  a  measured  error  and  in  all 
cases,  the  errors  will  be  substantially  reduced  so  the  system  would  quickly  acquire  a  zero 
wavefront  error  state  during  further  iterations.  It  is  important  that  the  system  reliably 
detects  the  null  state  of  the  wavefront  error  and  otherwise  maintains  its  proper  sign; 
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Fig.  7 

4.  Implementation  of  the  method  into  a  solar  wavefront  sensor 

A  wavefront  sensor  will  consist  of  two  masks  which  have  linear  independent,  preferr* 
ably  orthogonal,  shifts  A  and  A'.  A  pK>ssibility  to  generate  transparent  masks  is  given  by 
optical  modulators  commonly  used  in  photonics  research.  A  particularly  inexpensive 
approach  might  be  the  implementation  of  liquid-crystal  pocket  TV’s  (LCTV’s).  They 
make  good  optical  modulators,  with  20:1  contrast  in  white  light,  but  their  pixel  size 
might  be  somewhat  large  (0.1mm)  which  results  in  a  bulky  optical  setup. 

One  way  to  construct  a  LCTV-based  wavefront  sensor  is  shown  in  fig.  7.  A  polariz¬ 
ing  beam  splitter  is  used  to  generate  two  focal  planes  in  linear  polarized  light  for  the  two 
masks;  only  one  is  shown  in  the  figure.  A  modified  LCTV,  where  the  polaroids  have  been 
stripped  from  both  sides,  is  placed  in  the  focus.  A  piezo-elastic  modulator  (PEM)  is 
placed  before  the  LCTV;  the  modulator  serves  to  rotate  the  plane  of  polarization  by  90 
degrees  back  and  forth  at  a  high  rate,  say,  20  KHz.  The  liquid  crystal  rotates  the  plane  of 
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polarization  further,  depending  on  the  video  signal  at  the  pixels.  An  analyzer  is  placed 
following  the  LCTV  which  converts  the  polarization  modulation  into  modulation  of  inten¬ 
sity.  The  modulator  acts  as  a  contrast  inverter  in  conjunction  with  the  analyzer,  the 
video  picture  on  the  LCTV  is  constantly  contrast  inverted  at  the  modulation  frequency. 

Additional  optics  that  are  not  shown  reimage  the  telescope’s  entrance  pupil  onto  a 
diode  matrix.  Each  diode  samples  an  area  of  the  size  of  the  desired  subaperture,  there  are 
aa  many  discrete  diodes  as  there  are  subapertures  in  the  pupil.  Intensity  fluctuations  are 
measured  here;  because  of  the  modulation  the  intensity  fluctuations  are  offset  at  the 
modulation  frequency.  This  facilitates  their  detection  because  DC  offsets  in  the  liode 
preamps  are  thrown  out.  Lock-in  amplifiers  serve  to  demodulate  the  signals;  they  now 
represent  local  wavefront  tilts  and  they  are  sent  to  the  servo  system  which  controls  the 
active  device. 

The  LCTV  is  driven  by  a  video  camera  located  at  a  third  focal  plane  and  followed 
by  a  video  processor  which  generates  the  necessary  difference  pictures.  The  diff’  rence  pic¬ 
tures  are  generated  from  long  exposure  averages  which  are  updated  every  30  60  s. 

Alternatively,  the  analyzer  behind  the  LCTV  can  be  replaced  by  another  polarizing 
beam  splitter  followed  by  a  second  diode  matrix  in  the  pupil  image  in  order  to  save  light. 
This  approach  would  make  the  setup  somewhat  more  complex. 

5.  Signal  magnitude  and  noise  considerations 

We  shall  discuss  photon  noise  in  this  section  only.  The  required  signal  magnitude 
can  be  estimated  by  asking  how  many  photons  are  necessary  to  detect  and  maintain  the 
zero  wavefront  error  state  to  a  given  precision.  The  signal  for  wavefront  errors  close  to 
zero  is  given  by  eqn.  (7),  if  difference  measurements  are  made.  We  shall  now  make  esti¬ 
mates  for  the  various  factors  that  enter  this  equation  and  compare  the  signal  with  the 
bias  term  in  order  to  obtain  the  relative  intensity  variation. 

B:  average  mask  amplitude  transmittance  =  0.5. 

C:  mask  contrast;  depends  on  the  typical  rms  contrast  C^,  of  a  long  exposure  granula¬ 
tion  picture.  Under  average  seeing  conditions,  we  will  have  C,„,  =  0.08.  If  A  is  half 
the  correlation  scale  (0.7  arc  sec),  then  the  difference  image  will  also  have  about  3% 
rms  contrast,  or  10%  peak-to-peak.  Thus,  C  =  10  if  the  full  range  of  the  mask  is 
used. 
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I  F{al7)\ *:  power  integral  in  eqn.  (3).  The  integral  of  the  entire  power  spectrum  |  /'(*)|  * 
is  (O.OS)*  »  10~*;  a  factor  1/2  applies  because  the  integration  is  carried  out  over  half 
the  Fourier  plane  only.  The  sine  factor  in  eqn.  (3)  can  be  accounted  for  by  another 

factor  1/2  -*  |F(<r/2)|»  =  j  CiJ,,  «  2.5  x  10-‘. 

We  have  to  compare  the  result  P(i)  =  P(s)  -  P'(i)  with  the  average  pupil  intensity 

sum  I!  =  2fi*/.  Hence,  an  estimate  of  the  intensity  difference  is  given  with: 

=  2  X  10-»  (^(i)  -  ^(s  -  a/2)  )  =  2  X  10"*  A^(f)  (8) 

If  we  require  that  the  signal  is  photon  noise  limited,  we  obtain  a  minimum  number 

of  detected  photons  for  various  rms  wavefront  errors,  as  shown  in  table  I. 

We  now  make  an  estimate  of  the  number  of  available  photons.  The  area  of  a 
subaperture,  i.  e.,  the  area  subtended  by  a  single  detector  in  the  pupil  (cf.  fig.  7)  has  to  be 
at  least  so  large  that  the  target  structure  (granulation)  can  be  resolved.  Let  us  assume  a 
resolution  of  1.5  seconds  of  arc,  which  corresponds  to  a  subaperture  of  7cm  in  diameter. 
There  arc  1.7  x  10*  photons  A  i~‘  per  resolution  element  in  the  visible  at  about  SOOnm, 
this  number  increases  with  the  size  of  the  sensor  field  of  view.  If  we  allow  for  a 
wavelength  range  of  200  A,  sample  with  a  2KHz  rate  and  account  for  light  losses  in  the 
wave  front  sensor  (2  polarizers,  various  surfaces,  detector  quantum  efiiciency  etc.)  with  a 
factor  1/10,  we  arrive  at  total  photon  numbers  that  are  shown  in  table  II.  The  fraction  of 
light  needed  from  the  main  optical  beam  for  the  wavefront  sensor  working  photon  noise 
limited  is  shown  in  table  III. 


Tabic  I 


rms  WFE 

[radj 

P/S 

minimum  f  photons  (S/P)* 

A/10 

0.628 

l.SxlO”* 

5.9x10* 

A/20 

0.314 

6.3x10"* 

2.5x10* 

A/50 

0.126 

2.5x10"* 

1.6x10* 

Table  n 


sensor  field  size  f  photons  in  field  #  photons  per  measurement 


[arc  sec] 

[angitrom'^ 

200  A,  5x10-*  $ 

2x2 

5.0x10* 

5.0x10* 

5x5 

1.9xl0‘® 

1.9x10* 

10  X  10 

7.«xl0‘® 

7.6x10* 

Table  m 

sensor  field  size 

Fraction  of  light  used  for  WFS, 

at  rms  WFE 

[arc  sec] 

A/10 

A/20 

A/50 

2x2 

2x10-* 

8x10-* 

5x10-* 

5x5 

5x10-* 

1x10-* 

8x10"* 

10  X  10 

7xl0-‘ 

5x10-* 

2x10"* 

6.  Conclusions 

The  discussion  in  the  previous  section  shows  that  the  proposed  wavefront  sensor  is 
very  efficient.  Between  0.1%  and  5%  of  the  light  is  needed  for  the  sensor  to  operate 
photon-noise  limited  under  all  conceivable  situations.  Even  if  we  account  for  additional, 
instrumental  noise  the  amount  of  light  needed  still  might  be  about  as  low  as  10%. 

Some  comments  are  now  in  order  to  prevent  undue  criticism.  First,  effects  of  aniso- 
planatism  are  not  considered  here.  It  is  understood  that  anisoplanicity  poses  a  severe 
problem,  and  it  is  not  guaranteed  that  the  proposed  method  will  work  satisfactorily  if  the 
field  is  cut  down  to  approximately  the  correlation  scale  of  the  target  structure.  Much 
more  detailed  study  is  necessary  to  understand  all  the  effects  that  might  arise  from  such  a 
situation;  the  previous,  crude  analysis  is  certainly  insufficient.  But  the  technique  should 
work  with  fields  as  small  as  5  arcsec,  which  matches  the  typical  size  of  an  isoplanatic 
patch.  Similar  problems  will  arise  with  any  kind  of  wavefront  detection  using  solar 
granulation  as  target  structure. 

Second,  one  might  argue  that  a  subaperture  size  of  7cm  is  too  large  for  average  see¬ 
ing  conditions,  so  the  sensor  will  sucessfully  operate  under  good  seeing  conditions  only. 
One  has  to  keep  in  mind,  however,  that  wavefront  slope  measurement  are  made  and  pure 
wavefront  tilt  is  a  good  wavefront  error  approximation  to  a  7cm  region  if  r,  is  as  small  as 
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4cm!  (This  is,  by  the  way,  true  for  any  kind  of  slope  measurement).  In  fact,  the  pro¬ 
posed  technique  will  operate  over  a  large  range  of  seeing  conditions,  with  only  the  worst 
cases  excluded. 
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